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We consider the failure of localized control in a nonlinear spatially extended system caused by 
extremely small amounts of noise. It is shown that this failure occurs as a result of a nonlinear 
instability. Nonlinear instabilities can occur in systems described by linearly stable but strongly 
nonnormal evolution operators. In spatially extended systems the nonnormality manifests itself in 
two different but complementary ways: transient amplification and spectral focusing of disturbances. 
We show that temporal and spatial aspects of the nonnormality and the type of nonlinearity are 
all crucially important to understanding and describing the mechanism of nonlinear instability. 
Presented results are expected to apply equally to other physical systems where strong nonnormality 
is due to the presence of mean fiow rather than the action of control. 
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It has been known for a long time that systems de- 
scribed by nonnormal evolution operators (operators 
with non-orthogonal eigenfunctions) often display rather 
surprising dynamics. For instance, turbulence in shear 
flows often develops for Reynolds numbers where the ba- 
sic flow is still linearly stable. The critical Reynolds num- 
ber was found to depend rather sensitively on the geom- 
etry of the system and the roughness of the boundaries. 
Several studies [1-3] have linked the onset of turbulence 
to a nonlinear instability arising from the interaction be- 
tween the nonlinearity of the Navier-Stokes equation and 
the nonnormality of its linearization caused by significant 
mean flow. More recently the idea of a nonlinear insta- 
bility has been used to explain the disagreement between 
the predictions of the linear stability analysis and exper- 
imental data for the contact line instability in gravity 
driven spreading of thin liquid films [4]. Nonnormality 
can also arise in the absence of mean flow as a result 
of localized feedback control [5,6]. In this paper we use 
the idea of a nonlinear instability to explain the failure 
of localized control of a spatially extended nonlinear sys- 
tem in the presence of extremely weak noise. We extend 
and refine ideas described in [2,5] by incorporating the 
information about the spatial degrees of freedom. 

All studies of strongly nonnormal systems conducted 
up to now have analyzed the mechanism for nonlinear in- 
stability by concentrating only on the temporal dynamics. 
Although the importance of spatial degrees of freedom is 
generally recognized, the complexity of the problem usu- 
ally prevents consistent spatiotemporal analysis. As the 
subsequent discussion shows, spatial information is cru- 
cial for understanding the mechanism that leads to non- 
linear instability, which involves transient amplification 
of deviations produced by nonlinear terms. However, be- 
fore developing the spatiotemporal description, it will be 
useful to review some results of the temporal analysis. 
Following [2] let us consider a system 



where L is a stable linear operator and / is a nonlinear 
function of its argument. In a purely linear system the 
disturbance decays asymptotically in time. However, if 
L is nonnormal, this exponential decay can be preceded 
by a transient. The strength of nonnormality can be 
determined by the transient amplification factor 



\\m\ I 

7 = max - — -— -- — max 

' t,m 110(0)11 * I 



(2) 



such that 7=1 for L normal. The maximal transient 
amplification is achieved for the optimal initial distur- 
bance (/)(0) = 0opt at the optimal time t = topt [7]. 

Any initial disturbance (t>(0) with a nonvanishing com- 
ponent along (popt will be transiently amplified as well. 
For small disturbances, a = \\(j>{0)\\ ^ 1, the linear terms 
will dominate, so at the peak of the transient we will have 
||^</'(iopi)ll ~ ll?^(^opt)ll 7'^- For a quadratic nonlinear- 
ity; l!/(</'(*opt))ll (7^'')^; so it will produce an integrated 
deviation in (1) of order topti'^cr)'^ in the same amount 
of time it takes the linear operator to amplify the initial 
disturbance by 0(7). The temporal analysis makes an 
implicit assumption that generically this deviation has 
again a nonvanishing component along the optimal dis- 
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so it will be transiently amplified by L 
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in the same way as the initial disturbance. (As we will 
see below, this assumption can break down for spatially 
extended systems due to their high symmetry.) The devi- 
ation due to nonlinear terms will grow producing a posi- 
tive feedback loop, if iopt(7f)^ ^ cr, and decay otherwise. 
The critical magnitude of a disturbance needed to boot- 
strap the nonlinear instability will, therefore, scale like 
a ^ t'^ptl^^ for a quadratic nonlinearity. While in some 
systems like channel flow topt ~ 7, more often the depen- 
dence on topt is too weak to be of any importance, e.g., 
for both coupled map lattices [5,8] and partial differential 
equations [6] with localized control topt ~ log 7. In the 
latter case we have a simple power law scaling (7^7" 
with an exponent a — —2. 

Sometimes, temporal analysis is sufficient and does 
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FIG. 1. The noise level for which localized linear control 
fails as a function of the transient amplification factor for 
different types of nonlinearities (0^ - squares, 4>x4>- diamonds, 
- open triangles, 0^ - circles, df' - filled triangles). Straight 
lines are theoretical fits with slopes given by (5). 

give correct predictions for spatially extended systems. 
For instance, a controlled lattice of coupled quadratic 
maps [5] has indeed produced the scaling exponent a = 
—2. However, sometimes the predictions of the temporal 
analysis arc clearly wrong, suggesting that the spatial 
structure of disturbances plays an important role and 
cannot generally be ignored. A particularly simple ex- 
ample of the failure of temporal analysis is provided by 
the generalized Ginzburg-Landau equation (GGLE) 
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which (aside from the stochastic term ^) is of the same 
form as (1). The dynamics of GGLE can be made linearly 
stable via feedback control imposed at the boundaries 

(/-(0,t)=0, (t>'{L,t)= [ K{x)(j){x,t)dx, (4) 

Ja 

where K{x) is an appropriately chosen gain function. As 
we have shown previously [6], the application of spatially 
localized control (4) to the spatially extended system (3) 
makes the linearized dynamics strongly nonnormal. We 
therefore expect the nonlinear instability to play a promi- 
nent role in destabilization as large transient amplifica- 
tion makes the dynamics extremely susceptible to noise. 
Numerical simulations of (3) with a power law nonlinear- 
ity f{4i) oc and random noise ^ uniformly distributed 
on (—(7, ct) show (see Fig. 1) that the critical noise level 
resulting in the failure of linear control scales as a power 
law cr ~ 7", with an exponent well approximated by 



a = 



-n/(n- 1), 
-1, 



n = 2,4, 
n = 3, 5, 



(5) 



The exponent for n even is correctly predicted by a prop- 
erly generalized temporal analysis. Indeed, for controlled 
GGLE to'pt ~ log 7 ~ so comparing stochastic distur- 
bances of order a with the distortions of order (7(7)" pro- 
duced by a combination of transient growth and nonlin- 
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earity, we immediately obtain the scaling cr ~ 7 
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FIG. 2. The wavonumbers q„ and the corresponding eigen- 
values A„ = 1 — ?n for the original and controlled system. The 
system size here and in the rest of the paper is I = 20. 



However, the corresponding exponent is inconsistent with 
our numerical results for n odd. 

This discrepancy calls for the development of a more 
accurate theory that will be capable of explaining the ef- 
fects of arbitrary types of nonlinearities. In particular, we 
would like to know why an advective term (px^i produces 
the same scaling as a simple quadratic nonlinearity <jP 
despite their different symmetry properties. As it turns 
out, the explanation can be obtained rather easily by 
conducting a spatiotemporal analysis of the bootstrap- 
ping mechanism. Indeed, transient amplification repre- 
sents just one aspect of the nonnormal dynamics. The 
other aspect ignored by the temporal analysis is the fo- 
cusing of the initial disturbances in the direction of the 
most strongly nonnormal eigenfunctions. 

It is easy to sec that due to the translational invariancc 
of the linear operator L = 1 + d^, its eigenfunctions are 
sinusoidal, with or without control. The boundary condi- 
tions (4) determine the wavemimber q of an eigenfunction 
and the corresponding eigenvalue A = 1 — g^, such that 
the eigenfunction is stable when g > 1 and unstable oth- 
erwise. In particular, the eigenfunctions of the original 
system (no feedback, K{x) = 0) arc Uk{x) = sm{qkx) 
with qk = 7r(fc — l/2)/l, so at least one will be unstable 
for I > 7r/2. All eigenfunctions Vk{x) = sm{q'f,x) of the 
controlled system are stable with wavenumbers (/^ > 1. 
In general, might be complex, but we can always force 
them to be real. This is done throughout the paper by 
calculating K(x) using Linear-Quadratic control [9]. Ap- 
pearance of complex eigenvalues does not affect the fol- 
lowing analysis. As Fig. 2 shows, feedback (4) shifts all 
wavenumbers from the unstable band Q+ = (0, 1) into 
the stable band Q_ = (l,oo). The new wavenumbers 

cluster most tightly in a rather narrow band Qn cen- 
tered at q* ~ 1.36. As the size I of the system grows, 
an increasing number of eigenfunctions of the original 
system becomes unstable and gets squeezed into by 
feedback. As a result the distance between wavenum- 
bers in Qy shrinks and the corresponding eigenfunctions 
become increasingly aligned (nonnormal) [6]. 

Now, let us consider what happens with an initial dis- 
turbance (f){x,0). Let us first concentrate on the linear 
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FIG. 3. Fourier spectra of the first five adjoint eigen- 
functions of the controlled system. The expansion is in 
the basis {un} of cigenfunctions of the original system, 
Wkix) = Y,n ^kn am{qnx). 

effects. As the set {vk} is complete, in the absence of 
noise the general solution of the linearized equation (3) 
is given by 

oo 

= XI ^kVk{x) exp(Afcf), (6) 
fe=i 

where = 1 — ((/jj.)^. The coefficients Cfe can be found 
using the adjoint cigenfunctions Wk{x): 

{wk{x),Vk{x)) 

where we assume that the cigenfunctions are normalized 
such that ||w'A:|| = \\vk\\ = 1. As the eigenvalues are all 
stable, it is clear that transient amplification can only re- 
sult from large values of the coefficients Ck- For Ck to be- 
come large two conditions have to be satisfied. First, the 
numerator in (7) should not be small, i.e., the initial dis- 
turbance should not be orthogonal to the adjoint eigen- 
vector Wk- As adjoint eigenvectors satisfy the orthogo- 
nality condition {iUk,Vm) = for fc 7^ m, their Fourier 
spectra are localized to the unstable band Q+ (see Fig. 
3). Therefore, only disturbances with significant spectral 
content in (5+ will be transiently amplified. (Such distur- 
bances will grow indefinitely in the uncontrolled system; 
control makes this growth transient.) This is illustrated 
in Fig. 4 which shows the transient amplification for si- 
nusoidal initial disturbances: 

* mX,0)\\ ^(a:,0)=sin(gx) 

The second condition is that the denominator in (7) 
should be small, which can only happen for strongly non- 
normal cigenfunctions Vk with G Q||. As a result, 
the Fourier spectrum of the transiently amplified distur- 
bances will be strongly focused into the band This 
focusing effect can be clearly seen in Fig. 5 which shows 
the spectrum of the linearized GGLE driven by random 
noise. The spectrum is computed for the "worst case per- 
turbation," as this is the type of disturbance that leads 



FIG. 4. Transient amplification factor for sinusoidal initial 
disturbances as a function of their wavenumber. 

to both the failure of linear control and more generally to 
the onset of nonlinear instability. Specifically, the state 
is expanded in the basis {vk } 

00 

cj>{x,t) = Y,Ck{t)s\n{q'kX) (9) 
fe=i 

and the "worst case" spectrum Fk[(j)\ is obtained by find- 
ing the maximal values of |cfc(t)| for each k. The Fourier 
coefficients outside of Q\\ are seen to be exponentially 
small. These two intimately related aspects of the non- 
normal dynamics - transient amplification and focusing 
- are likely to be quite common in other strongly nonnor- 
mal spatially extended systems. For instance, a similar 
clustering of cigenfunctions is found in a model describing 
thermally driven spreading of liquid films [10]. 

Having understood the linear dynamics, let us now 
consider what happens when nonlinearities come into 
play. For sufficiently small a the nonlinearities will 
hardly change the linear dynamics. Their effect will 
be limited to "filtering" the transiently amplified distur- 
bances, changing their magnitude and spectral content. 
As the transiently amplified disturbances have a very nar- 
row spectrum, the spectrum of the signal produced by 
the nonlinear terms will also consist of several narrow 
peaks, as long as we consider nonlinearities of the power 
law type /((/>) = with moderate n. (High powers are 
not interesting as the scaling exponents (5) for the even 
and odd powers become indistinguishable. Besides, most 
physically interesting nonlinearities have low powers.) 

For instance, the spectrum of a quadratic nonlinear- 
ity, be it 0^ or (j)x<j), will only contain frequencies which 
are either sums or differences of frequencies g^, i.e., 0, 
\q'm~Q'n\^ "^^'rm ^^'^ ^'m+l'n- the Fouricr coefficient cor- 
responding to the frequency |g^±g^| is of order CkCm, the 
only significant contributions are produced when both q'^ 
and q'^ lie in Q\\. This means that the spectrum of the 
quadratic term will be localized near g = and q = 2q* 
(see Fig. 6a) . The disturbances with q « 2q* are strongly 
damped, so the primary effect of most any quadratic non- 
linearity will be to transfer the excitations from the band 
(5|| back into the band (5+, where they will again be tran- 
siently amplified. After one full cycle involving transient 
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FIG. 5. Normalized Fourier spectrum for a linear system 
driven by random noise uniformly distributed on {—a, a). 

amplification, focusing, and nonlinear filtering, an initial 
(e.g., stochastic) disturbance of order a will produce a de- 
viation of order (70-)^. Therefore, for (tct)^ > a the low 
wavenumber disturbances will be driven predominantly 
by the nonlinear term, bootstrapping a nonlinear insta- 
bility. A similar picture will be observed for /(0) = </>" 
with n = 4, 6, • • •. The spectrum of /(<?!>) will contain a 
strong component in and the bootstrapping will oc- 
cur for (70")'' ^ cr, so the critical noise will scale like 7", 
with a given by (5). 

The case of odd powers is substantially different. For 
instance, the spectrum of a cubic nonlinearity, f{(j)) — (p^, 
will only contain wavenumbers ± (7^ ± g^j. As the 
corresponding Fourier coefficients will be of order CmCnCk, 
the spectrum will be strongly localized near q = q* and 
q = 3q* . Furthermore, as Fig. 6 shows, the spectral peaks 
of the nonlinear terms do not broaden. Therefore, a cubic 
nonlinearity will not transfer excitations from Q|| back 
to and no bootstrapping will occur. The quintic 
nonlinearity, f{<p) = (j)^, is expected to produce similar 
results as its spectrum wiU be localized near q* , 3q*, and 
5^*, and so on. Destabilization will nevertheless occur 
for any power n when the nonlinear terms become of the 
same order of magnitude as the linear terms, (70")" ^ 70", 
i.e., when 7a = 0(1). At this point, the predictions of 
linear stability analysis become invalid. Therefore, for n 
odd the critical noise will scale like 7^^, justifying the 
second part of (5). The result for even powers is not 
changed, since nonlinear instability occurs for levels of 
noise much smaller than those at which linear stability 
analysis breaks down. 

Summing up, we can conclude that, at least for a 
simple equation such as the stochastically driven GGLE 
studied here, the failure of localized control can be ex- 
plained by a straightforward spectral analysis of transient 
dynamics. Moreover, spatiotemporal analysis appears to 
be crucial for understanding the mechanism of nonlinear 
instabilities in spatially extended systems in general. In 
particular, as the focusing effect described in this paper 
is an inherent feature of strongly nonnormal dynamics, 
its applicability is not constrained to the control prob- 
lem considered here. A similar analysis could provide 
valuable insights into stability and control of shear flows, 
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FIG. 6. Normalized Fourier spectrum of (a) the quadratic 
term, <f>'^, and (b) the cubic term, 0^. 



driven contact lines, and magnetic plasmas [11]. 
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